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Abstract 

Neutron stars provide a natural laboratory to test some unique implications of Quantum Chro- 
modynamics (QCD)- the underlying theory of strong interactions- at extreme conditions of very 
high baryon density. It has been suggested that the true ground state of QCD is strange quark 
matter, and, consequently, neutron stars should convert to strange quark stars under suitable 
conditions. Substantial efforts have been, and are being, spent in studying the details of such 
conversion. In this report, we show that the presence of high magnetic field, an essential feature 
of neutron stars, strongly inhibits the conversion of neutron stars to bare quark stars. 
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Neutron stars (NS) are compact objects formed as an aftermath of supernovae, resulting 
from the gravitational collapse of massive stars which have exhausted their nuclear fuel. The 
central density of NS is believed to be very high (6 — 10 times nuclear saturation density, 
depending on the total mass). They are also characterized by very high magnetic fields 
(10 10 — 10 12 G) at the surface. While the origin of such magnetic fields in NS remains 
unclear (trapping of all magnetic lines of force during the collapse [l| can barely yield 
10 8 — 10 9 G), the observed slow-down rate of pulsars, which are rotating NS, do require such 
magnetic fields. Thus, the canonical picture of the classical pulsar mechanism involves [l| 
a magnetic dipole at the centre of a rotating NS. The maximum mass of NS is somewhat 
model dependent, but in any case, of the order of the Chandrasekhar limit. 

In the classical picture, i.e., before the advent of QCD, it was believed that the only 
compact astrophysical objects were the white dwarfs and NS; anything other than these two 
classes would end up in black holes. However, QCD predicts the existence of a new class 
of compact stars, called quark stars [2] (QS), which again obey a maximum mass limit 
of the order of Chandrasekhar limit. It was proposed [4], Q some time ago that the stable 
strange quark matter (SQM) was the true ground state of QCD and as such, all NS ought 
to evolve, under suitable conditions, to this more stable configuration of QS, more correctly, 
strange stars (SS). Some tentative SS candidates are the compact objects associated with 
X-ray bursters GRO J1744 - 28 6], SAX J1808.4 - 3658 Q and X-ray pulsars Her X - 1 
8j. It has also been argued in the literature that sub- millisecond pulsars, if they actually 
exist, must be QS; the NS would not be stable at such high rotation speeds. 

There could of course be primodial SS, which could provide a natural explanation of the 
cosmological dark matter (CDM) within the standard model 9|. For the present purpose, we 
do not address them but rather concentrate on the evolution of NS to SS. There are several 
plausible scenarios where NS could convert to quark stars (ultimately SS), through a "seed" 
of external SQM [lOj], or triggered by the rise in the central density due to sudden spin-down 



It has been argued in the literature that such a conversion would not only provide a 



natural engine for the gamma ray bursts (GRB) but can also explain their observed beaming 



property 
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Starting from the original work of MIT group [10l . Il3j , a number of authors has studied 
the details of the conversion of NS to QS/SS; they have been summarized in our previous 
chapters and works 1J, ll5|]- (For the sake of brevity, it is not possible to cite all of them 
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here.) Suffice is to say, most of them have ignored some essential features of NS, namely the 
rotational effect, the general relativistic (GR) effect and, most importantly, the topic of this 
chapter, the role of the strong magnetic field. We have systematically tried to incorporate 
these effects in a series of works. In Ref. 14, we have shown that the conversion process is 
most likely a two-step process. The first step involves the conversion of NS to a two-flavour 
QS, occurring on millisecond timescales. The second step comprises the conversion from 
two-flavour quark matter to the stable three-flavour SQM, through weak interaction, the 
corresponding time scale being of the order of 100 seconds. In our more recent work [15j , we 
have exhibited the major role played by the GR effects in the rotating stars; the conversion 
fronts propagate with different velocities along different radial directions, a finding which 
could not have been anticipated from newtonian or special relativistic (SR) analyses. We 
saw that the velocity of the conversion front increased outward and engulfed all the matter, 
converting the NS to a bare QS. 

In the present work, we report the startling finding that incorporation of the magnetic 
field drastically alters the scenario. The formalism of the calculation, as reported in Ref. 15, 
remains unchanged and we need not discuss them in detail here; only the bare essentials are 
mentioned below. 

We start with the metric [3] 

ds 2 = -e~< +p dt 2 + e 2a {dr 2 + r 2 d8 2 ) + e^- p r 2 sin 2 e{d(f) - udt) 2 (1) 

describing the structure of the star, with the four gravitational potentials a,7, p and uj, 
which are functions of 9 and r only. The Einstein's equations for the potentials are solved 
through the 'rns' code with the input of an equation of state (EOS) and a central 
density. The nuclear matter phase is described by the non-linear Walecka model [18( and 
the quark phase by the MIT Bag model ljj. We, as usual, treat both the nuclear and 
quark phases as ideal fluids. The operative equations, then, are the continuity and Euler's 
equations; (for details, please see Ref. 15) 

1 (de de\ 1 / dv dv \ 2v v ^ (d^ da\ 

tu \dr dr J W 2 \dr dr J r r \dr dr ) 

I (dp dp\ 1 (dv dv\ 1 ( j^pl -£pP q^u) j\ 

w \dr dr J W 2 \dr dr J 2 y dr dr dr J ' 

where, v is the r.m.s. velocity, w is the enthalpy and W is the inverse of the Lorentz factor. 
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The presence of magnetic field alters these equations in a straightforward manner. The 
continuity equation is unaffected but the Euler's equation picks up an additional term, 



corresponding to the magnetic pressure ^ [2fJ. In the above Euler's equation, we replace 
pressure p by pi — p + With this modification, the calculation proceeds exactly as in 
Ref. 15. 

In keeping with the canonical picture, we assume that the magnetic field of the NS is due 
to a dipole at the centre of the star. At any point (r, 9), the field due to the dipole in the 
radial direction is given by 

where is the permeability in free space and m is the magnetic moment. The dipole is 
assumed to be along the polar direction and 9 is the angle with the vertical axis. We focus 
mainly on the range B sur f ace ~ 10 8 - 10 12 G. 



We now proceed exactly as in 15] . We define v as the front velocity in the nuclear matter 
rest frame and n = is the square of the effective sound speed in the medium, dr and dr 
are connected by the relation 



where G is given by 



_ ei + P - ei~Pr 2 uj 2 sin 2 9 , . 

G = y ^ ( 6 ) 



The other parameters in Eq. 2 and 3 are defined as 

vursin9 2u 2 rsin9 2uj 2 e i ~ p rsin9 

A = ~ cT~ " 1; E= Cl + Al ; 

Bl vursin9 2uoe 1 ~ p r 2 sin9 vrsin9 

Al = e 7+p - e J ~ p r 2 u 2 sin 2 9; 
Bl = -e 7+p - e^- p r 2 u 2 sin 2 9- 



where 



Cl = Vr 2 u> 2 sin 2 9 — e 2p \ 

After a bit of algebra, very similar to that in Ref. 15, we can write the generic differential 
equation for the velocity v of the conversion front 

dv W 2 v[K + Kl - K2] 



dr 2[v 2 (l + G) 2 - + v 2 G) 2 } 



(7) 
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FIG. 1: Variation of velocity of the conversion front along the equatorial direction of the star for 
three different cases, namely SR, GR and GR+M. 



where, 



K = 2n(l + v 2 G) [-L + — + - + 

\ or or r r 
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Integrating eqn. (5) over r from the center to the surface, we obtain the propagation velocity 
of the front as a function of r. 

We perform our calculation for a NS with a central density of 7 times the nuclear satura- 
tion density, (the corresponding Keplerian rotation rate being 0.89 x 10 4 sec _1 ). For the sake 
of comparison, we also study a slowly rotating star, with rotational velocity of 0.4 x 10 4 sec _1 . 

In fig. 1, we show the propagation velocity of the conversion front along the equatorial 
direction for the Keplerian star. The solid curve denotes the SR calculation, the dotted 
curve the GR result (both reported earlier 15J) and the dash-dotted curve the result of 
the present calculation, including the GR and magnetic field effects (GR+M). The surface 
magnetic field is ~ 5 x lO 10 ^. Fig. 2 shows the relative comparison between the slow star 
and the Keplerian star. Self consistency condition yields that the surface magnetic field for 
the slow star is slightly higher ~ 5.5 x 10 10 G. 
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FIG. 2: Variation of velocity of the conversion front along the equatorial direction of the star 
rotating with two different velocities for two cases. GR and GB.+M. 
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FIG. 3: Variation of velocity of the conversion front along the radial direction (polar and equatorial) 
of the star rotating with Keplerian velocity. 

From figs. 1 and 2, we clearly notice that near the centre of the star, the magnetic field 
produces a 'braking' effect. Only at large radial distances does the GR effect start to dom- 
inate and accelerate the conversion front. The effect of the rotational speed is quantitative; 
we find qualitatively similar behaviours both for slowly rotating and Keplerian stars. Fig. 3 
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FIG. 4: Variation of the velocity of the conversion front showing the 'braking' effect due to the 
magnetic field for different cases. 

shows the velocity of the conversion front along the polar and equatorial directions for the 
Keplerian star. x = cos9, defined w.r.t the vertical axis of the star. 

This 'braking' effect is naturally due to the additional pressure of the magnetic field. Our 
detailed calculation shows that for surface field less than 10 10 G, the 'braking' effect on the 
conversion front is not at all discernible; the GR effect dominates from the very beginning 
and the entire NS may get converted to a bare QS/SS. For higher magnetic fields, well within 
the canonical range, the effect however is much more drastic, as shown in fig. 4. 

Fig. 4 shows the results for both slowly rotating as well as Keplerian stars, along the 
equatorial radii. We see that for a surface magnetic field of 5.5 x 10 10 G along the equatorial 
direction, the velocity drops sharply and becomes almost zero at a distance of 3km (where 
equatorial radius is 15km) from the centre for the Keplerian star. Comparable and qualita- 
tively similar conclusions are also obtained for the slow star. In the polar direction, magnetic 
filed necessary to stall the front is slightly smaller than that in the equatorial direction. If 
the magnetic field is of the order of 10 12 G at the surface, the quark core shrinks within a 
distance of ~ 100m. For still larger fields, the front fails to start propagating. 

The implications of these results are indeed far-reaching. We are forced to conclude that 
the conversion to QS/SS is strongly inhibited in classical NS characterized by large magnetic 
fields. In particular, it is obvious that while these NS can at best convert to hybrid stars, 
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with a small quark core, the magnetars (NS with very high magnetic fields) have no chance 
to convert to QS/SS. Only very young NS, carrying small magnetic fields, have a chance to 
convert to QS/SS. 

For completeness, we would like to mention that the conclusion of this study does not 
depend much on the assumption of a dipole magnetic field. We have also employed an 
alternate field configuration which has been used by some authors 21[ in recent times. The 
conclusion remains unaltered. 
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